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Abstract: We consider several hints for new physics involving CP-asymmetries in 
B-decays and interpret them in terms of generic contributions to effective Wilson coef- 
ficients. The effects we focus on are: the differences in the fitted value of sin 2/3 versus 
the ones directly measured via the time dependent CP asymmetries in B — > J/ipK or 
via B — > ((f),r]')K; the difference between the direct CP asymmetries in B~ — > K~tt 
and B° — > K~tt + and the ~ 2.2<t indications for the CP-asymmetry in B s — > J/ip<f). To 
alleviate concerns regarding the disagreement between inclusive and exclusive V u b, we 
show that our results hold even without the inclusion of V u b in the analysis. We find 
that no matter what kind of new physics (NP) is invoked to explain these effects, its 
effective scale is bounded from above from a few hundred GeV to a few TeV depending 
on specific assumptions regarding the type of new physics. The only exception to this 
is when the NP contribution is assumed to reside entirely in LR operators in K mixing, 
then the scale of NP can be as high as around 24 TeV; however, this case cannot ac- 
count for CP asymmetry in B s — ► J/ip4> or a difference in sin 2(3 from penguin modes 
compared to that from J/ipK or for that matter the large difference seen between direct 
CP asymmetries in K~n + and in K~tt°. 



1. Introduction 



The only source of flavor changing interactions in the Standard Model (SM) is provided 
by the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix. Unitarity imposes testable 
constraints on the magnitude and phases of its elements; in particular, the relation 

v td v t t + v cd v; b + v ud v^ = o (i.i) 

has received considerable attention because of its strong sensitivity to the single CP 
violating phase that appears in the CKM. The extraction of magnitudes and phases of 
the various terms in Eq. (1.1) from a large number of observables is complicated by the 
presence of hadronic uncertainties. Due to the superb performance of the two asymmet- 
ric B-factories, it has been established that the Standard Model's CKM-paradigm [1] 
works to an accuracy of around 15-20% [2]; therefore the effects of New Physics (NP) 
are expected to be a perturbation and sub-dominant. The search for NP therefore re- 
quires very good control over theory errors and high statistics data from experiments. 
The combination of very precise experimental results from the B-factories BaBar and 
Belle and from the Tevatron experiments CDF and DO, with recent progress in lattice 
QCD, namely the improved determination of Bk [3] and of the B — > 7r form factor [4], 
leads to the emergence of several indications of possible deviations from the SM. We 
focus on the following issues: 

(a) The deduced value of sin 2(3^ differs from the directly measured value at the 2a 
level. The observables that are used to deduce the value of sin 2/3 are: the determina- 
tions of |T4b| and \V c b\ from inclusive and exclusive b — > (u,c)£i/ decays, the (indirect) 
CP violating parameter ek, and the ratio of the meson-antimeson mass differences in 
the B s and B d systems. In this deduction of sin 2(3 the only CP- violating quantity is 
that in the neutral kaon system. If the SM description of CP-violation through the 
CKM-paradigm with a single CP-odd phase is correct, then the value of sin 2(3 thus 
obtained should agree with the directly measured value of sin 2(3 in B-factory exper- 
iments. Following this logic, we will term this deduced value of sin 2(3 as sin2/9sM- 
Actually, since there is some ~ 2a discrepancy between the values of \V u b\ extracted 
from inclusive and inclusive decays [5], there is considerable motivation for obtaining 
sin2/3sM without using \V ub \ which we call sin2/3 s 1 ^ ub . Of course sin2/9g^ ub will only 
be of use if it has reasonably small errors. That this has become possible, due to the 
improved determination of Bk from the lattice, was recently empahsized in [6]. 

There are two important ways for extracting directly sin 2(3 via measurements of 
time dependent CP asymmetries. First there is the gold-plated {i.e. free from hadronic 

t\Ve recall the usuage of two equivalent notations: (<j>i , (fa , fe) = (P, ce, 7). 
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uncertainties to a very high degree of accuracy) measurement of sin 2(3 via the time de- 
pendent CP asymmetry in B — > (J/tpK) [7]; for clarity we will denote this as sin 2(3^ K . 
A second way to measure sin 2/3 is via b — > s penguin transitions such as 5 decays to 
t/Xs, i^s^s, n°K s , p°K s , uK s , f K s , n°n°K s , <jm°K s , K + K~K° [8,9], etc. 
Unfortunately, this method has some hadronic uncertainties. In the original papers 
this was crudely estimated at pa A 2 ~ 5% [9, 10]. In the past few years these modes 
received considerable theoretical attention [11-16] and as a result of that it now seems 
that amongst the two body modes, the r]'K s and <pK s final states receive hadronic cor- 
rections at the few percent level and are therefore very clean. For this reason, in the 
present work we will only include these two penguin modes and we will term sin 2(3 
extracted from their weighted average as sin 2(3^')k 

If the CKM description of CP violation is correct then all three determinations 
of sin 2(3 should agree with each other. In fact, both the "predicted" values of sin 2(3, 
whether one uses \V u b\ or not, differ from the directly measured values by ~ 2a. 

(b) New physics in b — > s penguin amplitudes is also hinted at by the fact that 
sin 2(3^k differs from sin 2(3u jT) i)k mentioned above by around 1.5a. 

In passing, we want to briefly mention that there is another feature of the time 
dependent CP-asymmetry measurements in various penguin modes''' that deserves dis- 
cussion. While the difference between Sj/^,k s and S^nguin, f° r eac h of the penguin 
modes is not that significant, a specially intriguing feature is that for many modes the 
central values of the asymmetry tends to be smaller than Sj/^k s - Since sin 2(3^k is less 
than sm2(3sM, this obviously implies that the central value of sin 2/3 pengu i n for almost 
all modes is also smaller than sin2f3 S M- 

It is also useful to recall that comparison of Sj/^k s with S^^k indicates whether 
or not b — > s penguin transitions are affected by NP. Since the time dependent CP in 
both tree and penguin modes necessarily involves B^B^ oscillations, comparison between 
sin 2(3^k and sin 2/5(0, rf) K does not tell us anything about whether there is new physics 
in Bd mixing. 

To the extent that sin 2(3^k differs from sin 2(3 sm, the possibilty of NP contributions 
to B d and/or K mixing (emphasized in particular in Ref. [17]) cannot be ruled out. 
Furthermore, once NP is invoked, we need to be careful in identifying which observables 
are sensitive to the type of NP that may be out there. 

(c) Another hint that b — > s penguin transitions may be exhibiting a non-standard 
CP-odd phase comes from the comparison of the partial rate asymmetry in B° — > 
K + tt~ and B + — > K + tt°. Experimentally this difference has been determined to be 
14.4±2.9% [5]. These two decays are closely related as they simply require switching the 

tTime dependent CP asymmetries in a generic B — ► / mode are denoted with the symbol Sf. 
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spectator (u,d) quarks. Therefore, the difference between these asymmetries vanishes in 
the limit of exact isospin and should be small. In sharp contrast, experimentally the two 
asymmetries are found to have an opposite sign and the result 14.4% is non-vanishing 
by over 4 a. It is difficult to rigorously assess the full significance of this unexpectedly 
large difference since we cannot reliably calculate, in a model independent fashion, the 
expectation from the SM taking QCD fully into account. In the QCD factorization 
approach [18,19], the predictions for BR and CP asymmetries in hadronic two-body 
B decays suffer from very large hadronic uncertainties, rendering problematic their use 
in NP searches. However, see for instance Refs. [20], a case can be made for NP in 
the difference AA CP = A CP (B- -> R-n ) - A CP (B° -> K~k + ). Ref. [20,21] show 
that in the QCDF analysis of this quantity, most parametric uncertainties that occur 
for individual asymmetries cancel out and the theoretical prediction becomes under 
reasonable control (see Ref. [22] for an alternate point of view) yielding (2.2 ± 2.4)% 
which is still about 3.5a away from the measured value. 

(d) The time dependent CP asymmetry in B s — > J/ip<p is free from hadronic 
uncertainties, has been recently measured at CDF [23] and DO [24, 25] and deviates 
from the SM at the 2a level [26,27]. Beyond the SM, a CP-odd phase in B s mixing is 
required in order to explain this discrepancy. 

(e) As is well known, over the past decade or so methods have been developed 
that allow direct measurement of all three angles of the unitarity triangle, (a, f3, 7) [2]. 
What makes these methods so atractive and useful is that attempts are made either not 
to use any theoretical input or assumptions or make the minimal use if necessary. For 
the angle f3 (via time dependent CP studies of e.g. B — > ipK s [7]) and the angle 7 (via 
direct CP studies in B^ — > (DK, D*K, DK*) and/or time dependent CP measurements 
in B°,B° — > (DK,D*K,DK*) [28-33]) no theoretical assumptions are needed. The 
resulting precision is largely data driven with an irreducible theory error of 0(0.1%) for 
(3 and < 0.1% for gamma [2]. For extracting the angle a a simple time dependent study 
of B°, B° — > 7r + 7r~ does not suffice and an isospin analysis [34-36] becomes necessary 
entailing a somewhat larger irreducible theory error, 0(few %). Currently, the angle a 
is being extracted by using branching ratios and CP asymmetries in B — > (im, pp, pn). 

Thus another useful avenue to exploit in order to test the CKM-paradigm and to 
constrain NP is to fit the Unitarity Triangle utilizing only the three angles (a, (3, 7) 
which are directly measured without theoretical assumptions or input. We will use this 
approach to extract the resulting values of the Wolfenstein parameters p and rj and 
compare them with the values deduced from the use of other methods and inputs. In 
particular, we show that once again the use of appropriate lattice matrix elements for 
ex, AM s ^, V c b and with or without V u b leads to ~ 2a deviations in one or both of the 
Wolfenstein parameters. 
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(f) Finally, we mention another approach to determine the fitted value of sin 2/3 
that does not use V u b determined from exclusive or inclusive semi-leptonic decays, rather 
this makes use of the directly measured values of a and 7 (see point (e) above) along 
with €k and the ratio of mass differences in B s and Bd mesons (see point (a) above). 
In the tables in Fig. 2 we list different ways of arriving at the fitted values of sin 2(3 
and indicate the resulting tensions. 

These anomalies in CP asymmetries involving B and B s mesons that we mention 
above are all at the (2 3) a level and may be indicative of new physics in Bd, B s 
and/or in K mixing and also in b — + s penguin transitions. Note that this new physics 
necessarily has to carry with it a beyond the SM CP-odd phase as all the observables 
being discussed here involve CP violation. 

In this paper we analyze these hints for new physics from an effective theory point 
of view: we parametrize NP contributions to various operators in terms of effective 
scales that, with the aid of these measurements, turn out to be constrained from above 
and from below. The interpretation of the two sigma discrepancies described above in 
terms of upper bounds on the scale of new physics has a caveat. In our approach NP 
contributions are essentially proportional to 1/A 2 where A is a generic NP scale. The 
tensions that we discuss translate into 1/A 2 > (and hence A < 00) at the 2a level; 
therefore, we are able to set an upper limit only at the 68% (and maybe 95%) C.L.. 
Beyond this confidence the measurements we are considering are compatible with the 
SM and no upper bound on NP is implied. 

The paper is organized as follows. In Sec. 2 we summarize the present state of 
the various fits that we use. In Sec 3 we perfom a model independent analysis of new 
physics in B d /K mixing and in b — > s penguin amplitudes. In Sees. 4 and 5 we interpret 
the discrepancies in the fit to the unitarity triangle (UT) and in penguin amplitudes 
in terms of NP contributions to various operators. In Sec. 6 we summarize and discuss 
our findings. 



2. Present status of the SM fits 

The set of inputs that we use in the fit is summarized in Table 1. at^^p,^ and the 
direct determinations of a that we obtain from the isospin analysis of B — > (titt, pp, pn) 
and B — > pp decays, respectively (we use the latter when discussing for NP effects in 
mixing so as to avoid pollution from possible NP contributions from b — > d penguins) . 
The direct determination of 7 is taken from the model independent UTfit analysis of 
B -> £>(*)#(*) decays [37,38]. The explicit expressions for AM Bq , AM B JAM Bd and 
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(41.67 ±0.68) x 10~ 3 incl [27] 
\V cb \ = { (38.7 ± 1.35) x 1(T 3 excl 
(41.0 ± 0.63) x 10~ 3 comb 

Am Bd = (0.507 ±0.005) ps" 1 [27] 
AA CP = (14.8 ± 2.8)% [27] 
<x™, mP x= (88.7 ±4.7)° 
77! = 1.51 ±0.24 [41] 
V2 = 0.5765 ± 0.0065 [43] 
?73 = 0.47 ± 0.04 [45] 
r) B = 0.551 ± 0.007 [46] 
S^k s = 0.672 ± 0.024 [27] 
S v , Ks = 0.59 ± 0.07 [27] 
£ = 1.211 ±0.045 [48] 
B K = 0.720 ± 0.013 ± 0.037 [3] 
A = 0.2255 ± 0.0007 [49] 



(39.6+^) x 10~ 4 incl [27] 
\V ub \ = { (33.8 ±3.5) x 10" 4 excl [4] 
(37.4 ±2.1) x 10- 3 comb 

Am Bs = (17.77 ± 0.10 ± 0.07) ps" 1 [40] 

7 = (78 ± 12)° [37,38] 

a pp = (87.8 ± 5.6)° 

m tiPole = (172.4 ± 1.2) GeV [42] 

m c {m c ) = (1.224 ± 0.057) GeV [44] 

e K = (2.232 ±0.007) x 10~ 3 

= ( 0.76 •!;;;;• V -2.37l°; 33 ) [25] 
S*k s = O.Ut^l [27] 



f Bs \JB Ba = (0.304 ±0.032) GeV [47] 

k £ = 0.92 ± 0.02 [17] 

f K = (155.5 ± 0.2 ± 0.8 ± 0.2) MeV [27] 



Table 1: Inputs used in the unitarity triangle fit. Quantities not explicitly given are taken 
from Ref. [5]. 



£k in the SM are: 



A Mu„ = 2 ' UU = lig^gll = ^m' w m Bq f Bq B Bq rj B S (x t ) \V tb V t tf , (2.1) 



m B a 



12tt 2 



AM, 



m Bd 



V t 



td 



(2.2) 



£k = 



G F m wfK m K £ 



B k k £ Tm(jnS (x e ) (V cs V c * d ) 2 + 2 m S (x c , x t )V cs V; d V ts V t * d 



+V2S (x t )(V ts V; d f) . (2.3) 
The B parameters (for M = K, B d , B s ) parametrize the matrix elements 

(M\QM\M) = lm 2 M f 2 M BWM (2.4) 
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Figure 1: Unitarity triangle fit in the SM (68% C.L.). The solid black, solid blue and dashed 
red contours are obtained using (a, (3, 7), (ek, AMb s , AMfl d , V c b) and (ek, AMb s , AMs d , 
V c b, V u b), respectively. 

where the operator Qi is given in Eq. (4.2). The Bm parameters are renormalization 
group invariant quantities and differ from the corresponding B^ s (p) by a perturbative 
factor (see, for instance, Refs. [46,50] for the details of this standard procedure). The 
quantity n e comes from the inclusion in Ek of the 1 = component of the K — > 7T7t 
amplitude [17,51-54]. The loop-functions can be found, for instance, in Ref. [17,50]. 
In this paper we will not concern ourselves with possible NP contributions to EW 
operators in the Kaon sector, whose effect is to alter the extraction of the factor n £ 
from data on e'/e (See Ref. [51] for a complete discussion of this issue). 

Our fitting procedure consists in writing a chi-squared that includes all experimen- 
tal measurements and lattice determinations. This procedure implies that all system- 
atic errors are treated as gaussian. While the true nature of systematic uncertainties 
remains subject of debate (see, for instance, the prescriptions adopted in Refs. [55,56]), 
we believe that our choice is preferable to flat systematic pdf's for several reasons: 
gaussian systematics lead to more conservative determinations of confidence level in- 
tervals; moreover, systematic errors in both lattice QCD and experiments are usually 
obtained by combining multiple sources of uncertainties, thus partially justifying our 
assumption. As usual we adopt the Wolfenstein parametrization of the CKM matrix 
and truncate the expansion at 0(A 4 ). The 68%C.L. allowed regions in the (p,fj) plane 
are shown in Fig. 1. In Fig 2 we summarize the numerical results that we obtain for p, 
fj and sin 2/3. In order to better illustrate the anatomy of the 2a tension we decided to 
present the results corresponding to the inclusion of different sets of observables in the 
fit (we show also the result of the complete fit for comparison) . It is interesting to note 
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that the model independent determinations of a and 7 affect the fit only in absence of 
the 1 14b I constraint. In the figure we also include a graphical representation of the dis- 
crepancy between direct and indirect (SM prediction) determinations of sin 2f3 as well 
as a pull table in which we quantify this discrepancy in terms of standard deviations. 
In the pull table, the w/out {with) V u f, column refers to the treatment of V u b on top of 
a fit that includes Ek, AM Bg , \V c b\, a and 7; the reference values of the SM predictions 
that we use in the pull table are therefore: 

J 0.846 ±0.069 w/out \V ub \ 
Sm(2/3) = 1 0.747 ±0.029 withal " (2 ' 5) 

These results summarize nicely the 2a tensions (a) and (b) that we discussed in the 
introduction. 



3. Model independent analysis 

The results of the previous analysis can be interpreted in the context of new physics 
contributions to 5^-mixing (M^ d ) , Ek and to b — > s penguin amplitudes (A^ s ). For 
the sake of simplicity we consider only the two extreme scenarios in which we admit new 
physics effects to {M± 2 d , Ab^, s ) and (ek, A^ s ), respectively. In this section we adopt 
very general parametrizations of possible new physics contributions; the connection to 
actual mass scales will be discussed in Sees. 4 and 5. 

3.1 Scenario I 

We assume all new physics effects to be effectively taken into account via the introduc- 
tions of two extra phases, (pa and 6 a- 

< = «)sm ^ d > (3- 1 ) 

•A-b^s = (^U^s)sM 6 A ■ (3-2) 

The expressions for the time-dependent CP asymmetries become: 

S^k = sin[2O3±0 d )] , (3.3) 
S {M)K = sin[2(/3±0 d ±0 A )] . (3.4) 

Note that we are implicitly assuming that NP effects in b — > s penguin amplitudes are 
identical in the (f> and 77' modes. This is necessary in order to use a simple parametriza- 
tion as in (3.2). This assumption will be relaxed in the operator level analysis presented 
in Sec. 5 where we adopt QCD factorization. Furthermore, the extraction of 7 from 
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P 


V 


sin 2/3 


e K ,AM Bd3 ,V cb 


0.236 ± 0.063 


0.478 ± 0.066 


0.885 ± 0.082 


e K ,AM Bds ,V cb ,V ub 


0.152 ±0.028 


0.383 ±0.019 


0.749 ± 0.030 


e K ,AM Bds , Vcb,a,j 


0.209 ±0.049 


0.442 ± 0.046 


0.846 ± 0.069 


e K , AM Bd s , Va, a, 7, V ub 


0.150 ±0.024 


0.382 ±0.018 


0.747 ±0.029 


a, 7, S^k 


0.116 ±0.027 


0.341 ±0.016 




e K , AM Bdsl V cb) V ub) a, 7, S$ K 


0.127 ±0.020 


0.357 ±0.012 





«, AM e? , |V cb | m 0.885±0.082 

es-, AM i? , |V cb |, a, y ^ 0.846±0.069 

AM V IVobl, a, y, I V ub | SM I ■ o.747±0.029 



b->ccs tree ■ 0.672±0.024 



0K» > I 0.44!™ 

ijK° penguin (clean) ^ o.59±o.07 

(0,V)K k-I I 0.57±0.065 



K S K S K S 1 . 1 0.74±0.17 

n a K° 1 — • — 1 0.57±0.17 

p K s 1 • — 1 0-o3_ 02 i 

<"Kj penguin (other) i • 1 o.45±o.24 

foKs 1 — Jh H 0.62^q ]3 

j*A s 1 . 1 -0.52±0.41 



mode 


w/out V u b 


with 


Si/jK S 


2.4 a 


2.0 


S^Ks 


2.2 a 


1.8 cr 


S V 'K S 


2.6 a 


2.1 a 


S(<t>+ V >)K S 


2.9 a 


2.5 cr 



07T K 5 



sin(2fi) 



Figure 2: Results of the fit to the unitarity triangle within the SM. In the table on top we 
collect the results we obtain for different selection of inputs. The plot is a graphical comparison 
between the SM predictions given above and the direct determinations in b — ► ccs and b — > s 
penguin modes. In bottom-right table we show the deviation of the clean sin 2(5 measurements 
from the SM predictions obtained using ek, AM Bq , \V cbl a, 7. The last column shows the 
impact of \V ub \. 



B — > D^K^ decays is controlled by tree-level decays and is assumed to be insensitive 
to new physics effects. This assumption does not hold for a. In this case the isospin 
analysis extracts and effective angle given by a e fr = en — <pd + ^penguin- Here (pd is the 
same angle that appears in Eq. (3.1), while # pe ngum is a possible new physics phase in 
b — > d penguin amplitudes. In order to simplify the analysis we will utilize only the 
B — > pp channels in the extraction of cv e ff because, in this case, the penguin contribution 
turns out to be experimentally very small and we are justified in setting #p en g U i n = 0. 
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Figure 3: Model independent analysis of new physics effects in Bd (left panel), K (right 
panel) mixing and A\,^ s (both panels). Solid and dashed lines corresponds to the la and 2a 
contours, respectively. 



In this scenario, the theory prediction for sin(2/3) is obtained by excluding from the 
chi-squared both S^k and a; we obtain: 

f 0.867 ±0.080 without | Kb | ,„ K , 

Sm(2/?) = \ 0.747 ±0.029 with 1^1 " (3 ' 5) 

In the left panel of Fig. 3 the contours define regions whose projections on the axes 
yield the one dimensional ranges at 68% (solid) and 95% (dashed) confidence level. 
Without the inclusion of V u \, we obtain: 

<j> d = - (8.7t|£) (3.6) 
e A = -(3.8 ±2.1)° (3.7) 

As can be seen from Fig. 3 the negative error on <pd is highly non-gaussian. The 
sharp cut-off at <pd ~ —11° is due to the interplay between ex and AM^/AM^ on 
one side and S^k and B —>■ pp on the other (we remind the reader that the former 
extract /3 eff = (3 + (pa and the latter a cS = a — 4>d)- The inclusion V u b into the fit 
lowers considerably the predicted value of sin(2/3) (as can be deduced from Eq. (2.5) 
and Fig. 1), thus impacting strongly the extraction of <pd- On the other hand, 9a is 
essentially determined by the difference between the time dependent CP asymmetries 
in B — ► J/i^K on one side and B — > (<f),r] f )K on the other; hence the outcome of the 
fit, Eq. (3.7), is largely independent on the inclusion of V u \,. Numerically we find that 
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both the (f>d central value and error decrease by a factor of two: 

fa = - (3.2±1:1)° (3.8) 
9 A = -(3.5l?;jj)° (3.9) 

In this scenario we interpret the tension in the fit to the unitarity triangle in terms 
of NP contributions to the time dependent CP asymmetries in b — > ccs and b — > sss 
modes. The discrepancy between the predicted and "measured" values of sin (2/3) is 
explained by new physics contributions to Ba mixing; the difference between S^k and 
S(<j>,r)')K is induced by a new phase in the b — > s penguin amplitude. 

3.2 Scenario II 

We now assume the absence of new physics contributions to B& mixing and investigate 
the possibility that the tension in the fit is induced by new effects in K mixing. The 
discrepancy between the time dependent CP asymmetries in the b — > ccs and b — > 
s penguin modes still requires an independent NP phase. We adopt the following 
parametrization: 

e K = C £ (e K ) SM > ( 3 - 10 ) 

•A S ss = ("4-sss)sm e A ' (3-11) 

Note that the expressions for the time-dependent CP asymmetries become: 

S c5s = sin [2(3] , (3.12) 
S sSs = sm[2((3 + e A )\ . (3.13) 

In the right panel of Fig. 3 we show the 68% C.L. (solid) and 95% C.L. (dashed) allowed 
region in the (C e , 6a) plane. The one-sigma ranges for these two parameters without 
the inclusion of V u \> in the fit read: 

C £ = (1.31 ±0.14) , (3.14) 
9 A = -(3.6 ±2.3)° . (3.15) 

The impact of V u b shifts only slightly these values: 

C £ = (1.28 ±0.13) , (3.16) 
9 A = -(4.1 ±2.3)° . (3.17) 

In this scenario, the sin(2/3) prediction coincides essentially with S^ K and does not 
depend much on the inclusion of V u b, hence the amount of new physics required to 
bring Ek in agreement with the rest of the fit is quite insensitive to the V u b constraint. 



- 10 - 



4. Operator analysis of new physics in the fit to the UT 



From our previous discussion it is clear that the tension in the fits of the unitarity 
triangle are related to the presence of new physics either in Bd or K mixing. The 
effective Hamiltonian that describes meson mixing (B d , B s and K) can be written as: 



Glml 



n cS = {v tq v*,y [ J^CiOi + ) + h.c, (4.i) 



J=l i=l 



where we have (q, q') = (bd) , (bs) , (sd) for Bd, B s and K mixing. The operators are 
defined as follows: 

01 = (QlI^l) (QlI^l) (4.2) 

2 = (q' R q L ) (M (4.3) 
Os = (<i%q() (fe) (4.4) 
4 = (q' R q L ) (q' L q R ) (4.5) 
O, = fegf ) (qlqi) (4.6) 



and 0^2,3 are obtained from 01,2,3 via the L <-> R substitution. In the following we will 
consider new physics contributions to the Wilson Coefficients of the operators Oi and O4 
(O4 has the same anomalous dimension and matrix element as O4, hence we are really 
constraining O4+O4). Oi is the only operator that receives a non-negligible contribution 
in the SM. Contributions to O4 are especially interesting because they are enhanced by 
QCD running effects and by chiral factors [57] that appear in the calculation of their 
matrix elements; in particular, in the K mixing case these effects result in a two order 
of magnitude enhancement. We parametrize new physics contributions to the various 
Wilson coefficients as: 

**V> - -c^rv ■ (4-7) 

where we retained a factor l/(167r 2 ) to take into account a possible loop-suppression 
of NP effects* and we decided to factor out the CKM couplings. The factor —1 is 
introduced because we know from the model independent analysis of Sec. 3.1 that 
the required NP phase has to be negative. Combining Eqs. (4.1) and (4.7), the NP 
contribution to the effective Hamiltonian is 

(V t V*,) 2 e i{ c 

SK„ = - [ -^L ij O, + h.c. (4.8) 



"This factor appears in the denition of the effective Hamiltonian (4.1). 
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and the scale A absorbs every NP coupling, mass scale and loop function apart from 
the CKM matrix and the typical l/(167r 2 ). Using Eq. (4.7) we obtain Cnp/Csm — 
— e*^(700 GeV/A) 2 ; hence for A ~ 700 GeV the new physics and SM contributions to 
the Wilson coefficients are of similar size. In the remainder of this section we show the 
bounds on A that we obtain in the two scenarios we introduced in Sec. 3. We will also 
consider the possibility of simultaneous NP contributions to (Bd, B s ) and (Bd, B s , K) 
mixing 

4.1 New physics in Bd mixing 

We begin by assuming that new physics contributes to a single operator relevant to Bd 
mixing (this corresponds to Scenario I of Sec. 3.1) and take 5 NP Cf d ^ 0. We do not con- 
sider NP contributions to the other (bd)(bd) operators that appear in Eq. (4.1) because 
their hadronic matrix elements are all very similar (e.g. no large chiral enhancement 
of LR operators) and the outcome of the analysis does not change appreciably. Note 
that we introduce new physics to the Bd mixing amplitude only; as a consequence, the 
phase <fid is non vanishing and, since we do not allow contributions to B s mixing, the 
ratio AMbJ AM Bd will be affected too. Using the expression Eq. (4.7) for the generic 
NP contribution to C? d we have: 



M- 



12 



K 2 | SM e 



2i/3 



/«A< 



(4.9) 



<p(°) 



<5Cf rf =- 



Gp myj A 2 



without V, 



with V, 




where n = G^m^S^Xt). From the model 
independent results we obtained in Sec. 3, 
we see that the phase <pd = l/2arg(l — 
e^/nA 2 ) is negative, implying ip > 0. In 
Fig. 4 we show the allowed regions in the 
A, ip plane. We write a chi-squared that 
contain all the observables we discussed in 
Sees 2 and 3 (namely e K , AM Bda , V cb , \V ub \, 
7 from DWif" decays, a from B — ► pp, 
and S^jk s ) an d minimize with respect to all 
variables (including p and fj) . The contours 

are such that their projections on the axis correspond to the one-dimensional 68% CL. 
regions for A and (p. The green (dashed) and blue (solid) contours are obtained with 
and without the inclusion of \V u b\ in the fit, respectively. The presence of the upper 
limit A < 2.3 TeV reflects the two sigma effects Eqs. (3.6) and (3.8). The lower bound 



500 1000 1500 2000 25c 

A (GeV) 

Figure 4: New physics in Bd mixing. 



- 12 - 



A > 1. TeV is a direct consequence of NP contributions to X s a'. 



AM Bs ( AMr 



1 



l 



AM Bd \AM B J SM • (410) 

The qualitative impact of the \4& constraint can be inferred from the analysis of Sec. 3.1 
and from Fig. 3. A reduction in the absolute size of <pd translates into smaller values 
for if — ir. Finally, it is interesting to extract the predicted value for sin (2(3): 

sin(2/ 3 ) = ( - 82±01 ° * ithootV - . (4.11) 
V 1 ! \ 0.73 ±0.03 with V ub V ; 

The comparison of this result with the SM prediction given in Eq. (2.5), shows that 
in this scenario the tension between sin(2/3) and the CP asymmetries in the (pK and 
rj'K channels is somewhat eased: S^ +V /) K deviates from sin(2/3) at the 2.2/2.1 a level 
with/without the inclusion of V ub . 

4.2 New physics in both and B s mixing 

In this section we modify the analysis of Sec. 4.1 by allowing simultaneous identical 
NP contributions to Bd and B s mixing: 5Ci d = 5C*f s . This approach is inspired by 
a Minimal Flavor Violating (MFV) ansatz* in which NP contributions to are 
identical up to CKM factors: 

5M d l2 ^-{V tb V; d f e ^ 2 and 5M{ 2 oc - (V tb V:f ^ . (4.12) 

It is important to stress that we introduce this complex contribution only in flavor 
changing operators involving quarks. In explicit NP models that implement this idea, 
e.g. extra Z' with flavor changing quark couplings, it is important to keep under control 
phases that appear in interactions that conserve flavor or involve leptons (the latter, 
in particular, are constrained by CP asymmetries in exclusive b — > s£ + £~ decays). The 
inclusion of a complex correction to M{ 2 is of great interest because it allows to reconcile 
the constraint coming from the time dependent CP asymmetry in the B — > J / ip<p system 
with the rest of the fit. In the numerics we utilize the HFAG combination of DO and 
CDF data: <p s = — 0.76±g;|| or S = -2.37^;p [27] (information from the flavor specific 
life time and the semileptonic asymmetry is included). 

In this scenario the mass differences AM Bs , AM Bd , as well as the B q mixing phases 
are affected but the ratio AM B J AM Bd does not receive contributions. This is because 
the relative impact of NP on the B q systems is identical. The measurements that 
we include into the fit are Ek, AM Bds , V cb , \V u b\, 7 from D^K^ decays, a from 
B — > pp, S^k s an d Sj/^. The result of the analysis is summarized in Fig. 5. A 

^However, we stress that if the B-CP anomalies we discussed in here are confirmed, that would be 
quite inconsistent with the general notions and expectations of models based on MFV. 



-13- 



striking feature of these results is the large impact that the inclusion of V u b has on 
the allowed regions in the (A, ip) plane. The impact of V u b is to require a smaller 
4>d phase because the discrepancy between the direct and indirect determinations of 
sin(2/5) decreases in absolute value. On the other hand the observed discrepancy in the 
phase measured in the J/ip 4> system still points to large effects. The friction between 
these two competing effects results is responsible for the large shifts in the contours 
obtained with and without V u b- The predictions that we obtain for sin(2/?) are: 



sin(2/3) 



0.88 ± 0.07 without V ub 
0.75 ±0.03 with Kb 



5C, 



A 1 



In this case, the tension between the global extraction of sin(2/3) and the 
asymmetries is unaffected. 



4.3 New physics in K mixing 

The scenario described in Sec. 3.2 corre- 
sponds to new physics contributions to the 
K mixing amplitude only. We implement 
this framework by allowing contributions 
to either Of or Of : all new physics effects 
are confined to Ek- Note that this time we 
consider separately possible NP effects in 
the LR operator Of: because of the large 
QCD running effects on its Wilson coeffi- 
cient and of the chiral enhancement of its 
matrix element, the bounds that we ex- 
tract for this case are about one order of 
magnitude stronger then the ones we ob- 
tain for new physics in 0\. The explicit 



(4.13) 
rf)K CP 



<p( u ) 




without V u b 

with V ul 



A (GeV) 

Figure 5: New physics contributions to B d 
and B s mixing. 



formula for Ek that we use to study NP contributions to the Wilson coefficients of the 
operators Of and Of is obtained from the SM formula (2.3) via the substitution: 



S (x t ) 
S (x t ) 



So(x t ) 
S (x t ) 



o'V 



AvA 2 

B 4 K 44 



, for NP in C{ 



K 



3m 2 K 



kA 2 B k t] 2 4 ( m s + m d f 



for NP in Of 



(4.14) 
(4.15) 



= x 



where k = G 2 F m 2 w S (x t ), m s = (100 ± 10) MeV [58,59], m d = (4.6 ± 3) MeV [58], 
_E>4 = 1.03 ±0.06 [60] and K 44 ~ 3.7 is calculated below. The factor x, whose numerical 
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estimate is given below in Eq. (4.19), quantifies the different impact that new physics 
contributions to the Wilson coefficients C\ and C 4 have on e K and is enhanced by chiral 
and QCD-running effects. The masses of the strange and down quarks are defined in the 
MS scheme at the scale \xl ~ 2 GeV and are taken from Ref. [58,59] (the actual value 
of m s that we adopt reflects the dispersion of several lattice results). B 4 = B^ s (^i L ) 
is the bag parameter of the operator O4 and it has been calculated in quenched lattice 
QCD; the value we use is taken from Ref. [60]. Note that we have: 

(Of (/xl)> _ 3 m\ B 4 

(Of(/i L )) 4(m s + m d )2 B K ' 1 ' 1 

Finally, QCD effects in the running of the Wilson coefficients between fi H and \xl are 
summarized in the matrix K rs : 

^ = E(^ S) +^V 1 ' (4 - 17) 

i 

aw = ^^ s c s w, (4.18) 

s 

where 77 = a s (M^)/o; s (m t ) and the magic numbers a«, 6, and q have been calculated 
for — 2 GeV in Ref. [60]. In this analysis we take fi H = m t . The MS-bar scheme 
dependence of the bag parameters Bi can be removed by including part of the QCD 
running effects into their definition, thus leading to the introduction of the hat param- 
eters Bi. Eq. (4.15) is justified because we have B K K U = B K r] 2 . Putting everything 
together we find 

The presence of NP in M^ 2 affects the extraction of sin(2/3) from S^k- In the present 
case NP contributions to K mixing are proportional to (V t dV t * s ) 2 ; hence their effect is 
0(1) on ek ~ ImM/l but only 0(0.1%) on ReM^. In the following we will neglect 
such corrections to S^ K . In Fig. 6 we show the results of the analysis. The projections 
of the blue (solid) and green (dashed) contours onto the A and ip axes correspond to 
the one-dimensional 68% C.L. ranges. 

Note that the inclusion of V u b in the fit does not impact appreciably the allowed 
regions in the (A, ip) plane. In this scenario the role of NP effects is to reconcile the 
Ek constraint with the rest of the fit (V^, V u b, AM Bq , a, 7 and, especially, S^k)- 111 
particular, the inclusion of S^k in the fit renders the latter quite insensitive to V u b- In 
general V u b tends to slightly improve the overall consistency of the fit to the UT within 
the SM, therefore after the inclusion of Kb the bounds on A become slightly weaker. 
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Another interesting feature of Fig. 6 is the absence of a lower limit on the scale 
of NP. This happens because Ek is given by the imaginary part of M/|; hence a given 
correction can be obtained for any arbitrarily small A by appropriately choosing (p. 
Note that the NP contribution to Mj^ is proportional to — e^ 2/3+ ^ / A 2 and in the limit 
ip — > {it — 2(3, 2n — 2(3) the correction induced on ek vanishes: this feature is evident in 
Fig. 6 in which the asymptotic values of ip in the limit A — > are very close to 7r — 2/3 
and 2tt - 2(3. 

The upper bounds on the NP scale that we extract are about 1.8 -j- 1.9 TeV and 
23 -T- 24 TeV for the Cf and Cf scenarios, respectively. 



5. Operator analysis of new physics in b — > s amplitudes 

In this section we interpret the difference between the time dependent CP asymme- 
tries S^k and Sfarf in terms of new physics contributions to the QCD or EW penguin 
operators. The effective Hamiltonian responsible for the B — > ((p,r] f )Ks amplitudes is: 

An ( 6 6 \ 

V \ i=l i=3 / 

The definition of the various operators can be found, for instance, in Ref. [61]. Here we 
focus on two operators whose matching conditions are are likely to receive new physics 
contributions: 

Q 4 = (sL^T a b L ) £ (q lfl T a q) . (5.2) 
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Figure 7: New physics in b — > s penguin amplitudes. 



Q 3Q = (w»b L )Y,Q q (<rtv<i) • (5-3) 



We adopt the following parametrization of new physics effects: 



6a 



a* 



4,3Q 



A 2 4:71 [y/2 



(5.4) 



where we kept the coupling and loop suppression typical of QCD and EW penguins 
(the factor a Sje /(47r) = (e 2 , g 2 )/(167r 2 )). Note that we have absorbed all new physics 
couplings in the effective scale A. It is important to notice that once we introduce 
new physics contributions in these penguin operators we induce important effects in 
the B — > Kit system as well; in particular, we will consider the difference of CP 
asymmetries AAcp = Aqp{B~ — > K~tt°) — Acp{B° — > K~tt + ). In order to describe 
the impact of the new physics coefficients onto the <f>K, rj'K and Ktt system, we follow 
the QCD factorization analysis of Refs. [11,18,19]. Within the SM we find: 

S3^ K = 0.03 ± 0.01 , (5.5) 
SS^k = 0.01 ± 0.025 , (5.6) 
AA^p = (2.2 ± 2.4)% . (5.7) 

The errors in Eqs. (5.5-5.7) are obtained by varying simultaneously all the hadronic 
inputs that we take from Refs. [11,18,19]. For what concerns the B — > Kir asymmetries 
we quote also the results that we obtain for the various topological amplitudes: C/T = 
(0.16 ± 0.08) - z(0.08 ± 0.05), P E w/T = -(0.67 ± 0.03) - i(0.004 ± 0.06), \T/P\ = 
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0.24 ± 0.03. T and C correspond to the color allowed and color suppressed matrix 
elements of the tree- level operators Oi )2 ; P and P E w are the one loop matrix elements 
of the QCD (0 3 _ 6 ) and EW (O 7 _i ) penguin operators, respectively. In our conventions, 
the definitions of these amplitudes include the magnitude of the corresponding CKM 
factors but not their phases. 

In the left and right panels of Fig. 7, we present the allowed regions of the (A, tp) 
plane in presence of new physics contributions to C 4 and C 3 q. The blue and red shaded 
regions are obtained using the constraints from Su^k an d AAqp, respectively. The 
black streets cixc obtained by requiring both constraints simultaneously. The excluded 
regions within the AAcp contours correspond to a part of the parameter space that 
yields a too large value of AA C p. The irregular behavior of the S^^) K contours is 
due to the complicated structure of the theoretical errors on these quantities. For each 
point in the (A, ip) plane we determine the theoretical error by varying all the hadronic 
inputs; the resulting two-dimensional error function is then utilized in the chi-squared 
fit. 

The most important result of this analysis is the existence of an upper limit on 
the effective scale of about 400 GeV (200 GeV) for new physics contributions to QCD 
(EW) penguin operators. 

6. Summary 

We discussed several anomalies involving CP asymmetries in B and B s decays. The 
measured CP asymmetry in B — > J/ipKs when compared with the SM prediction from 
the fits of the UT seems to be too small by about 15% and this hints to new physics in 
either or K mixing. The non-vanishing differences between the time dependent CP 
asymmetries in B — > J/ipKs and B — > ((f),r)')K modes monitors the presence of new 
physics in b — > s transitions. The latter is also hinted at by the direct CP asymmetries 
in the Kir system (AA C p). The large asymmetry in B s — > J/ifxfi also indicates a non- 
vanishing beyond the SM phase in B s mixing on the face of a negligible asymmetry in 
the SM. 

The tension in the fit to the UT can, for example, be explained with an extra phase 
in Mf 2 whose value is found to be <pd ~ — (3± 1.5)° (we obtain — (9±3)° if no use of V u b 
is made) or by new physics in Ek for which we find C £ ~ 1.3 ± 0.1, where in the SM, 
C e — 1. The anomaly in the ((f), rf) system points to a new phase in b — > s amplitudes 
for which we obtain 9a ~ —(4 ± 2)°. 

These results can be interpreted in an effective Hamiltonian formalism in terms of 
NP contributions to Wilson coefficients. In this way we can translate these hints for 
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NP into scales at which we expect to find accelerator signals. Our results for different 
physics scenarios are summarized in the following table: 



Scenario 


Operator 


A (TeV) 


¥>(°) 


Bd mixing 




\ 


f 1.1 -5-2.1 no Kb 
k 1.4 -5- 2.3 with Kb 


\ 


f 15 4 92 no Kb 
_ 6 4 60 with Kb 


B d = B s mixing 


o[ d) & o[ s) 


\ 


f 1.0 -T- 1.4 no Kb 
k 1.1 -5- 2.0 with Kb 


\ 


' 25 4 73 no Kb 
_ 9 4 60 with Kb 


K mixing 


of 
of) 


< 1.9 

< 24 


130 -r- 320 






.25 -T- .43 
.09 4 .2 


4 70 
4 30 



Our main finding is that, no matter what kind of new physics is invoked to explain these 
effects, its effective scale is bounded from above at few TeV. The only exception are NP 
contributions solely confined to the LR operator in i^-mixing; however, it should be 
stressed that if NP affects only K mixing, then it cannot explain the difference in the 
extraction of sin2[3 from B — > J/ipKs and B — > (cf),r)')K and, in addition, it cannot 
account for both AAcp(kn) and the asymmetry in B s — > J/ipcp. 

Finally, let us comment on a similar analysis for the scale of NP presented in 
Ref. [62]. One important difference concerns the treatment of Ek'- we utilize the recent 
(2+l)-flavors determination of Bk from the RBC collaboration and the estimation of 
the effect of the 0-isospin K — > amplitude on K mixing (the factor k £ ) [17,51-54]. 
The combined effect of the updated values for these parameters is to strengthen the 
impact of the Ek constraint on the fit to the UT and to introduce the > 2a discrepancy 
responsible for the upper limits on the NP scales that we find. Another important 
difference of the present analysis from Ref. [62] is in the treatment of NP effects on the 
CP asymmetries in b — > s penguin and Kn modes, both of which are included in our 
analysis and not in Ref. [62]. 

Lastly, we want to briefly comment on the possibility of resolving these anomalies 
within the SM. In particular the impact of indirect CP violation in the K system, ex, 
depends crucially on the hadronic matrix elements Bk and on the precise value of |Kb| 
(we remind the reader that the p and 77 dependent part of Ek is proportional to |Kb| 4 )- 
While a rather large shift in a single input parameter (such as B K or Kb) is needed 
to reduce the discrepancy between the fitted and measured (via B — > J/^K S ) values 
of sin 2/5, a correlated set of smallish shifts in several inputs, while implausible, can 
certainly not be ruled out. However, for the effects that we discuss to disappear, any 
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such problems in the hadronic matrix elements from the lattice and/or elsewhere will 
only suffice, if simultaneously it is proven that the B s — > J/ip(p asymmetry and the 
smaller values of sin 2f3 from penguin modes were all a statistical fluctuation. 

Finally, let us summarize the impact of future experimental and theoretical progress 
on the anomalies we considered in this analysis. The most promising developments will 
be the high precision measurement of the B s — > J/ip<f) asymmetry, additional lattice- 
QCD calculations of Bk, the inclusion of 0(a^) and 0{a s /mi } ) corrections in the global 
fit to inclusive b — > civ decays for the extraction of \V c b\, and the calculation of the 
parameter k £ introduced in (2.3) using 2+1 flavors lattice-QCD (k £ is controlled by the 
matrix element of the QCD penguin operator Oq - not presently calculable with good 
accuracy with lattice QCD methods [63] - and can be extracted, within the SM, from 
the measurement of e'/e and the lattice determination of the matrix element of the 
electro- weak penguin operator 8 .). Because of recent progress in lattice calculations, 
the errors on Bk and k £ do not impact too strongly the fit to the unitarity triangle 
anymore; therefore new determinations of these parameters will serve mainly to build 
confidence in the central values that we are presently using. On the other hand, the role 
of \V c b\ is of the utmost importance: if the inclusion of higher order QCD corrections 
to inclusive semileptonic B decays will help closing the 2a gap between inclusive and 
exclusive determinations of |V^,|, the discrepancies we considered in this work will be 
strongly reinforced. Improved determinations of CP asymmetries in B — > Kit and 
b — > s penguin modes will most probably have to wait for LHC-b and/or super-B 
factories [64,65]; unfortunately the interpretation of these discrepancies relies heavily on 
QCD-factorization methods and suffers from our lack of control over power corrections. 
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